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ABSTRACT 

Recent analysis of Goddard High Resolution Spectrograph (GHRS) echelle spectra 
suggests ~ 30% variations in the D/H abundance ratio along the line-of-sight to 
the nearby (69 pc) hot, white dwarf (WD) G191-B2B (Vidal-Madjar et al. 1998, 
hereafter VM98). Variations in the D/H ratio on such short length scales imply 
non-uniform production/destruction of deuterium and an inefficient mixing of gas 
in the local interstellar medium (LISM). We reinvestigate the question of spatial 
variation of the local D/H abundance, using both archival GHRS spectra, and new 
echelle spectra of G191-B2B obtained with the Space Telescope Imaging Spectrograph 
(STIS) aboard HST. The STIS spectra were obtained in the high-resolution (E140H) 
mode and cover the wavelength region 1140 to 1700 A. Our analysis uses stratified 
line-blanketed non-LTE model atmosphere calculations to determine the shape of the 
intrinsic WD Lyman-a profile and estimate the WD photospheric contamination of the 
interstellar lines. Although three velocity components were reported previously towards 
G191-B2B, we detect only two velocity components. The first component is at v^ei ~ 
8.6 km s _1 and the second at v^i ~ 19.3 km s _1 , which we identify with the Local 
Interstellar Cloud (LIC). From the STIS data we derive D/H = 1.60±o.2? x 10" 5 for the 
LIC component, and D/H > 1.26xl0 -5 for the 8.6 km s _1 component (uncertainties 
denote 2a or 95% confidence limits). The derived D/H values in both components are 
consistent with (D/R) LIC = 1.5±0.1xl0~ 5 determined by Linsky (1998). The STIS 

2 NASA/Goddard Space Flight Center, Code 681, Greenbelt, MD 20771 

3 National Optical Astronomy Observatories, 950 North Cherry Avenue, Tucson, AZ 87519-4933 
4 Raytheon ITSS, NASA/Goddard Space Flight Center, Code 681, Greenbelt, MD 20771 
5 IACS/Department of Physics, Catholic University of America, Washington DC 20064 
6 Advanced Computer Concepts, Inc. /Goddard Space Flight Center Code 681, Greenbelt, MD 20771 
7 Dept. of Physics & Astronomy, University of Leicester, Leicester, LEI RH, UK 

8 University of Arizona, Lunar and Planetary Lab., West Gould-Simpson Bldg. Rm. 901, Tucson, AZ 85721 

1 Based on observations done with the NASA/ESA Hubble Space Telescope obtained at the Space Telescope Science 
Institute, which is operated by the Association of Universities for Research in Astronomy (AURA), Inc. under NASA 
contract NAS5-26555 



- 2 - 



data provide no evidence for local or component-to-component variation in the D/H 
ratio. 

Our re-analysis of the GHRS data gives essentially the same results as VM98, 
despite using two velocity components for the profile fitting (versus three by VM98) 
and using a more physically realistic WD Lyman-ct profile for G191-B2B. The GHRS 
data indicate a component-to-component variation as well as a variation of the D/H 
ratio in the LISM, neither of which are supported by the newer STIS data. The D I 
absorption in the GHRS spectrum is shallower than in the STIS spectrum. We believe 
the most probable cause for this difference in the two data sets is the characterization of 
the background due to scattered light in the GHRS and STIS spectrographs. The D/H 
ratios derived are sensitive to the background subtraction procedures employed. The 
two-dimensional MAMA detectors of STIS measure both the spatial and wavelength 
dependences of scattered light, allowing more accurate scattered light corrections than 
was possible with GHRS. 

Subject headings: white dwarfs (G191-B2B) — ISM: abundances — ultraviolet:ISM 

1. Introduction 

The best post-IUE measurement of the D/H ratio in the LISM within ~100 pc has been the 
HST-GHRS study by Linsky et al. (L93 & L95), who found D/H = 1.6 ± 0.2 x 10" 5 towards 
Capella, a nearby late-type spectroscopic binary system [d = 12.5 pc; (7,6) = (162°. 6, +4°. 6)]. 
Subsequent GHRS measurements towards other nearby late-type stars and WDs, reported by 
various authors indicate that within measurement uncertainties, the D/H ratio in the LISM is 
constant and the value is consistent with Linsky et al.'s results (Landsman et al. 1996, Piskunov 
et al. 1997, Dring et al. 1998). Only one result suggests significant variation of the D/H ratio 
within 69 pc. Using GHRS echelle data, VM98 reported the presence of three velocity components 
towards the white dwarf G191-B2B, with D/H ratios ranging between 0.9xl0" 5 and 1.56xl(T 5 
implying a variation of the D/H ratio by ~ 30% in the LISM. Towards the more distant star 
<50ri A {d ~ 350 pc), however, Interstellar Medium Absorption Profile Spectrograph (IMAPS) 
observations by Jenkins et al. (1999), suggest spatial variations of the D/H ratio. Current 
models of Galactic chemical evolution predict variations in the D/H ratio over length scales of 
~ 1 kpc (e.g. Tosi, 1998) but not over length scales as short as 69 pc. Possible implications 
of the spatial variation of the D/H ratio over such short length scales include a non-primordial 
source of deuterium production (Mullan h Linsky 1998). We have reinvestigated the question of 
spatial variation of the D/H ratio towards G191-B2B, using newer HST-Space Telescope Imaging 
Spectrograph (STIS) data. Compared to the GHRS echelle data, the STIS data have better 
scattered light corrections and include lines such as Si n (AA1190, 1193, 1260 and 1526) and Fe II 
A1608.5 which were not observed with GHRS. In this Letter, we report the D/H ratio derived from 
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STIS data and the disagreement between the D/H ratios derived from STIS and GHRS data . The 
most probable reason for this disagreement is the inadequate scattered light corrections available 
for the GHRS data (§6). An extended analysis of the STIS and GHRS data will be presented 
elsewhere. 

2. Observations and data reduction 

G191-B2B was observed with STIS on 1998 December 17 using the high-resolution E140H 
(R~110,000) mode in the ~ 1140 to 1700A region. These observations were part of STIS flux 
calibrations and the entrance aperture used was 0.2 x 0.2 arcsec. The use of this aperture 
decreases the effective spectral resolution, particularly in the bluer wavelength regions near the 
Lyman-a region (where the telescope PSF halo is more pronounced). The spectra were processed 
with the IDL-based CALSTIS reduction package developed by the STIS Instrument Development 
Team (IDT) at Goddard Space Flight Center (Lindler 1999). There was appreciable scattered light 
in the Lyman-a region (~ 11% of the continuum at 1213A) which was corrected using an iterative 
correction algorithm developed by the STIS IDT (Bowers & Lindler 1999, in preparation). This 
algorithm models several sources of stray light including the telescope PSF halo, echelle grating 
scatter and the detector halo. Archival GHRS echelle data of the G191-B2B (observed 1995 July 
26 to 28) were reduced using the IDL-based CALHRS routine developed by the GHRS IDT at 
Goddard. Three individual GHRS FP-SPLIT exposures are combined to obtain the final spectrum 
containing the D I and H I interstellar lines. The core of the saturated Lyman-a absorption dips 
below zero flux and was corrected by setting the inter-order coefficients a and b at 0.9 (Cardelli 
et al. 1993). Near the D I feature, the S/N per data point in the continuum for both GHRS and 
STIS is ~ 20, although GHRS has better wavelength sampling (0.003A versus 0.005A for STIS). 

3. Use of NLTE stellar atmosphere models 

To disentangle the interstellar D I and H I absorption from the observed profiles, it is essential 
to use the most physically realistic model for the intrinsic stellar Lyman-a profile. G191-B2B [d 
= 68.8 pc; (l,b) = (155°. 9, +7°.l)] (Vauclair et al. 1997), belongs to a class of hot, DA WDs that 
contain significant amounts of heavy elements such as C, N, O, Si, Fe and Ni in their atmospheres. 
Lanz et al. (1996) performed NLTE calculations including the effects of line-blanketing from 
more than 9 x 10 6 atomic transitions (mainly Fe and Ni) and matched the flux level and shape 
of the EUV spectrum of G191-B2B for the first time. The apparent effective temperature of 
WDs like G191-B2B is sensitive to assumptions about the photospheric composition (Barstow 
et al. 1998) and must be taken into account in modeling the stellar Lyman-a profile. Barstow 
and co-workers (in preparation) have refined their stratified line-blanketed NLTE calculations 
(Barstow et al. 1999) and the best-fit model atmosphere [T e jj = 54,000 ± 2000 K and log g = 
7.5± 0.03] is adopted in our analysis to predict the intrinsic WD Lyman-a profile and to check 
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for contamination of the interstellar lines by narrow WD absorption lines. The radial velocity of 
G191-B2B used in our analysis, estimated from STIS data of other WD lines (Bruhweiler et al. 
1999, in preparation) is 24.6 ± 0.4 km s _1 (including gravitational redshift). 

4. The number of velocity components in the line-of-sight 

In addition to the interstellar D I and H I absorption lines, the STIS echelle spectra show 
interstellar absorption due to N I (AA1 199.5, 1200.2 and 1200.7), C n A1334.5, C* n A1335.7, 
O I A1302, Si II (AA1190, 1193, 1260, 1304 and 1526), Si m A1206.5, Al n A1670.8, S n A1259.5 
and Fe n A1608.5. The interstellar N I A1200.7, Si n AA 1193 k 1304 and Fe n lines are not 
contaminated by WD lines. Figure 1 (a, b) shows the profile fits to the N I A 1200.7 and Si n 
A 1304 lines and two distinct velocity components are seen in these uncontaminated lines (all 
velocities are in the heliocentric frame). One component is at ~ 8.6 km s" 1 (hereafter referred to 
as comp 1). The other component is at 19.3 km s _1 , which is within measurement uncertainties 
of the projected velocity of the LIC (20.3 km s _1 ) in the line-of-sight to G191-B2B (Lallement et 
al. 1995). This component is also seen in the Capella data (L93), suggesting both the G191-B2B 
and Capella sightlines intercept the LIC. Figure lc shows the STIS and GHRS Si in 1206.5 A line 
utilized by VM98 to determine the number of components in the line-of-sight towards G191-B2B. 
Our GHRS profile is shifted to lower velocities by ~ 4 km s" 1 compared to the profile in the 
top panel of Figure 5 in VM98. Unlike VM98, we are able to obtain an excellent fit to both the 
STIS and GHRS profiles of Si in using only these two components. Detailed profile fitting of the 
other interstellar species confirmed that no more than two components are required (within the 
constraints imposed by S/N and spectral resolution of the STIS data) to yield acceptable fits. Our 
analysis of the STIS and GHRS data explicitly assumes the existence of two distinct components. 

5. Profile fitting of the interstellar D I and H I lines 

Each component is assumed to be homogeneous and characterized by a column density N, 
radial velocity v and a line-of-sight velocity dispersion defined by b = (2kT/m + ^ 2 ) 1 / 2 where ^ is 
the turbulent velocity parameter along the line-of-sight, T is the kinetic temperature and m is the 
ion mass. The D I and H I interstellar lines were fit simultaneously since they are separated by 
only 0.33A and the D 1 absorption is located on the wing of the broad H I absorption. Line profiles 
were convolved with either the STIS instrumental LSF for the 0.2 x 0.2 arcsec slit given by the 
STIS Handbook (Sahu 1999) or the two-component Gaussian LSF for GHRS given by Spitzer & 
Fitzpatrick (1993). The turbulent velocity parameters for the two components were determined 
by plotting the b values for the various atomic species as a function of ion mass m and performing 
a least-squares fit. The best-fit £ value for the LIC component is 1.7 km s^ 1 (consistent with 
L93, L95) while for comp 1, £ is 2.5 km s -1 . The STIS spectra near the Fe 11 (the heaviest ion) 
absorption have low S/N and the £ value is probably not very accurate. However, the derived 
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column densities are insensitive to the assumed values of £, and the D/H ratios presented here 
are not affected. For modeling of the Lyman-a profile, the velocities of the two components are 
kept fixed at 8.6 (comp 1) and 19.3 km s _1 (LIC) and the £ values are fixed at 1.7 (LIC) and 
2.5 km s _1 (comp 1). Two types of profile fits are done for the STIS and GHRS data sets: (1) 
keeping the value of D/H free in both components (STIS-FREE, GHRS-FREE) and (2) forcing 
the same value of D/H in both components (STIS-FIXED and GHRS-FIXED). Table 1 lists 
the results of the profile fitting for the STIS and GHRS data sets. The parameters obtained 
for the two components are the H I column density, N(H i), the temperature derived from the 
H I thermal velocity dispersion and the D/H ratio (listed for comp 1 and the LIC component 
respectively in columns 4 through 9). The uncertainties quoted for the D/H ratios denote the la 
(95%) confidence limits obtained using the method of constant x 2 boundaries (Press et al. 1992). 
Figure 2 (a, b) shows the best-fit models to the STIS and GHRS data for the STIS-FREE and 
GHRS-FREE cases respectively. The total H I column density towards G191-B2B obtained from 
the STIS data is ~2.04xl0 18 cm -2 , consistent with the value of 2.05xl0 18 cm" 2 derived from 
the best-fit parameters to the EUVE data over the wavelength rang e 100 to 500A (Barstow et al. 
1999). The D/H ratios derived for the two components from the STIS and GHRS data clearly 
disagree (compare STIS-FREE and GHRS-FREE cases in Table 1). 



6. Why do STIS and GHRS data give different values of the D/H ratios? 

Figure 3 compares the STIS and GHRS spectra in the region of the D I absorption. Note 
that the D I absorption in the GHRS spectrum is shallower than in the STIS spectrum. The 
difference in the derived D/H values is unlikely to be due to statistical fluctuations. For example, 
a value of D/H = 1.35 xlO -5 in both components is 3a above the best-fit determination from the 
GHRS-FIXED fit, and 3a below the best-fit determination from the STIS-FIXED fit. Another 
unlikely possibility for the difference in two data sets is time variability in the observed profile 
near the D I feature, perhaps due to a stellar wind. While Barstow et al. (1999) do suggest the 
presence of a weak stellar wind in G191-B2B to maintain the stratification of the Fe abundances, 
they point out that the wind must be less than 10 -16 M Q /yr to avoid elimination of the heavy 
elements in the photosphere. Such a weak wind would not be detectable, even in Lyman-a. The 
most probable cause is a systematic error in one of the data sets and two lines of evidence suggest 
that this error is more likely to be in the GHRS data. First, whereas x 2 l v value of 1.065 for the 
model fit to the STIS data indicates a 16% probability that the model is correct and that the 
uncertainties are correctly estimated, the x 2 l v value of 1.200 for the GHRS data indicates only a 
0.01% probability of this being true. Second, when the D/H ratio is kept fixed in both components 
(STIS-FIXED and GHRS-FIXED), the D/H ratios derived with the STIS data (1.71±o.mx10 -5 ) 
are consistent with the value of 1.5±0.1xl0 -5 determined for the LIC by Linsky (1998). The 
corresponding value derived for the GHRS data (1.17±o.'n xl0~ 5 ) is not consistent with observed 
LIC values. Fixed-pattern noise or wavelength drifts during FP-SPLIT subexposures in the 
GHRS data set could result in a shallower D I absorption. However, the observed scatter in 
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the flux level among the 144 subexposures (comprising the three FP-SPLIT GHRS exposures), 
shows good agreement with the errors predicted by the CALHRS routine, making this possibility 
unlikely. We believe the most probable cause for the difference in the two data sets around the D I 
absorption is the characterization of the background due to scattered light in the GHRS and STIS 
spectrographs. The D I line located on the wing of the H I profile, almost reaches zero flux in the 
core and is sensitive to the background subtraction procedures employed. The two-dimensional 
MAMA detectors and smaller pixel sizes on STIS allow a better estimate of the wavelength and 
spatial dependence of scattered light than the one-dimensional GHRS digicon science diodes. 
After the scattered light correction was applied to the STIS data, the residual flux in the core of 
the saturated H I profile is less than 1% of the continuum flux at 1213 A. The GHRS observations 
were done with the default STEP-PATT option, where the background is measured with the 
science diodes for only 6% of the total exposure time (Soderblom 1995). Due to the low S/N of 
this background measurement, only a low-order polynomial fit can be made to the background 
variations (Cardelli et al. 1990). Use of the standard values of the four echelle scatter correction 
coefficients recommended by Cardelli et al. (1993) yields a significant (~ 8% of the continuum 
flux at 1213 A) oversubtraction in the core of the saturated H I profile (see Fig 1, VM98). We 
corrected for this by multiplying the background by 0.9 prior to subtraction (§2), but if there are 
significant variations in the background with wavelength, this ad hoc correction is inadequate. 
Most previous D/H ratios derived with GHRS have used late- type stars as the background sources 
and these emission-line sources should have a much less significant scattered light problem. Given 
the importance of the question of spatial variations of the D/H ratio in the LISM, future STIS 
observations of G191-B2B with a smaller aperture (to achieve higher wavelength resolution), 
would be valuable. 

This research was supported by a GTO grant to the STIS IDT. MAB acknowledges support 
from PPARC, UK. We are grateful to Jeff Linsky, the referee, for insightful and helpful comments. 
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Table 1. RESULTS OF PROFILE FITTING 



Fit type 






N(H i) b 


Comp 1 

C 


D/H d 


N(H i) b 


LIC Component 

C 


STIS-FREE 1 


1.065 


457 


0.32±0.11 


11110±820 


2.24(> 1.26) t 


1.72±0.10 


6130±830 


STIS-FIXED 2 


1.065 


458 


0.41±0.04 


10570±520 


1 T-i i 0.32 
1.71±o.24 


1.65±0.05 


5740±800 


GHRS-FREE 1 


1.200 


740 


0.66±0.08 


8590±630 


0.66±o.ol 


1.51±0.08 


7090±780 


GHRS-FIXED 2 


1.203 


741 


0.53±0.04 


9040±510 


i.i7±8:i? 


1.62±0.04 


6900±770 



D/H d 



i.6o±8:i? 



1.29±8."i7 



a v = number of degrees of freedom 
b in units of 10 18 cm" 2 

c Temperature in K derived from H I thermal velocity dispersion 
d D/H in units of 10 -5 ; errors quoted denote 2a (95%) confidence limits 
Profile fitting keeping the value of D/H free in both velocity components 
2 Profile fitting forcing the same value of D/H in both velocity components 

tin the STIS-FREE profile fit, the region of the D I absorption line that contributes most to the D I 
column density of comp 1 reaches close to zero flux (see Fig 2a), resulting in a D/H ratio that is not well 
constrained. 
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Fig. 1. — Normalized STIS echelle spectra of interstellar absorption lines in the line-of-sight to 
G191-B2B are plotted against heliocentric velocity for (a) N I A1200.7 (b) Si n A1304 and (c) Si III 
A1206.5. The data points are shown with ± la error bars. The best-fit models to the absorption line 
profiles are plotted as continuous lines. The dotted lines indicate the normalized fluxes predicted by 
the WD NLTE stellar atmosphere models discussed in § 3. The lower panel for each absorption line 
shows the residuals to the fits with the y-axis in each lower panel ranging from -0.1 to 0.1 in relative 
intensity units. Both N I and Si n absorption lines are uncontaminated by WD photospheric lines 
and show two distinct velocity components (indicated by arrows), one at v^ei ~ 8.6 km s~ x (comp 
1) and another at Vh e i ~ 19.3 km s _1 (LIC component). In (c), the Si in GHRS profile is plotted 
as a histogram (± la error bars plotted every ~ 7 data points) over the STIS profile and shows 
the good agreement between the two datasets (§ 4). We are able to obtain excellent fits to both 
the STIS and GHRS profiles of Si III 1206.5 A using only two velocity components. 

Fig. 2. — Profile fits keeping the value of the D/H ratio free in both velocity components, obtained 
using (a) STIS data and (b) GHRS data. The observed data are shown as histograms. ± la error 
bars are shown for every ~ 25 data points. The best-fit models to the absorption line profiles are 
plotted as continuous thin lines and the lower panel for each absorption line shows the residuals 
to the fits. The individual component fit for the LIC absorption component is shown by dotted 
lines while comp 1 is shown by dashed lines The dot-dashed line shows the intrinsic WD Lyman-a 
profile predicted by the line-blanketed stratified NLTE atmosphere models (§3) which was used in 
the profile fitting. In the case of the GHRS spectra (b), the region ~ 1215.7 A which contains the 
geocoronal feature, has been removed from the data before profile fitting. 

Fig. 3. — A comparison of the STIS (continuous histogram) and GHRS (dotted histogram) data 
m the 1215 to 1216.5 A region. A 3-pixel smoothing has been applied to both data sets, and 
the zero levels and scale factors of both data sets have been adjusted to provide agreement with 
theoretical models. Note the D I line in the GHRS spectrum is shallower in comparison to the STIS 
spectrum. The D/H ratios for the two velocity components obtained from STIS data are consistent 
with (D/H)lic = 1.5±0.1xl0~ 5 obtained by Linsky (1998) and provide no evidence for a variation 
in the D/H ratio within ~ 70 pc. 
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